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ABSTRACT 

We present results from the mid-infrared spectral mapping of Stephan's Quintet using the Spitzer 
Space Telescope 10 . A 1000 km s -1 collision [t co i = 5 x 10 6 yr) has produced a group-wide shock and 
for the first time the large-scale distribution of warm molecular hydrogen emission is revealed, as 
well as its close association with known shock structures. In the main shock region alone we find 
5.0 xl0 8 M Q of warm H 2 spread over ^480kpc 2 and additionally report the discovery of a second 
major shock-excited H 2 feature, likely a remnant of previous tidal interactions. This brings the total 
H 2 line luminosity of the group in excess of 10 2 ergs -1 . In the main shock, the H 2 line luminosity 
exceeds, by a factor of three, the X-ray luminosity from the hot shocked gas, confirming that the 
H 2 -cooling pathway dominates over the X-ray. [Sin]34.82/im emission, detected at a luminosity of 
l/10th of that of the H 2 , appears to trace the group- wide shock closely and in addition, we detect 
weak [Fe n]25.99/im emission from the most X-ray luminous part of the shock. Comparison with shock 
models reveals that this emission is consistent with regions of fast shocks (100 < V s < 300 km s -1 ) 
experiencing depletion of iron and silicon onto dust grains. 

Star formation in the shock (as traced via ionic lines, PAH and dust emission) appears in the intruder 
galaxy, but most strikingly at either end of the radio shock. The shock ridge itself shows little star 
formation, consistent with a model in which the tremendous H 2 power is driven by turbulent energy 
transfer from motions in a post-shocked layer which suppresses star formation. The significance of 
the molecular hydrogen lines over other measured sources of cooling in fast galaxy-scale shocks may 
have crucial implications for the cooling of gas in the assembly of the first galaxies. 
Subject headings: galaxies: evolution - galaxies: individual (NGC 7318b; NGC7319) - galaxies: in- 
teractions - galaxies: intergalactic medium - physical data and processes: shock 
waves 



1. INTRODUCTION 

Stephan's Quintet (hereafter SQ) is a strongly inter- 
acting compact group which has p roduced a highly dis- 
turbed intragroup medium (IGM) (|Xu et al.ll2003l 120051) 
through a complex sequence of interactions and har - 
rassmcnt (jAllen fc SullivanlHoM IWilliams et ai]l200"2 ). 
This interplay has produced a large-scale intergalactic 
shock-wave, first observed as a narrow filament in the 
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radio continuum (jAllen fc 'Hartsuik erl Il972f) . and sub- 
sequently detected in the X-ray dPietsch et al.1 119971 : 
iTrinchieri etld][2005l IQ'Sullivan et al.N2009r i. The high- 
velocity (^1000 km s" 1 ) collision of an intruder galaxy, 
NGC 7318b, with the intergalactic m edium of the 
group ()Sulentic et alJl200lHXu et al.ll2003h is believed to 
be responsible for the shock-heating of the X-ray emit- 
ting gas. Optical emission line ratios and observed broad 
linewidths provide evidence that the region is powered 
by strong shocks and not star formation (Xu et al. 2003; 
Due et al. 2010, in preparation). 

The main elements of Stephan's Quintet are shown in 
Figure 1. Central to the system is the primary shock, 
as defined by the 20 cm radio continuum emission. NGC 
7318b, the intruder galaxy, lies t o the West o f the shock, 
while the large Seyfert 2 galaxy (jDurret]|1994[ ). lies to the 
East. Other members of the group are also indicated. 
The p eculiar extranuclear s tar formation region, named 
SQ-A (|Xu et al.lll999Ll2003f h lies at the extreme northern 
end of the main shock wave. NGC 7318a (west of NGC 
7318b) is also a strongly interacting group member, with 
NGC 7317 further away from the core. 

The unexpected discovery of extremely powerful, pure- 
rotational H 2 line em ission from the center of the shock 
( ApDlcton et al. 2006) , using the Spitzer Space Telescope 
(|Werner et al.l l2004L has sparked intense interest in the 
SQ system. Warm molecular hydrogen emission was 
found with a line luminosity exceeding the X-ray lu- 
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Fig. 1. — Schem atic composite im age of Stephan's Quintet consisting of a Galex Near-UV image (blue), an optical R-band 
image (cyan) from lXu et alj |l999), an IRAC 8/im image (yellow) and an image from the 1 6/xm IRS Peak- up Array (red). The 
contours represent the 20 cm emission associated with the shock as obtained by the VLA l|Xu et al.| [2003). The boundaries of 
the Spitzer IRS Spectral Mapping coverage (this paper) are shown as boxes. A color version of this image is available in the 
online publication. 



minosity from the shock. The mid-infrared (MIR) H2 
linewidth was resolved (a = 870 km s -1 ) suggesting that 
the H2-emitting clouds carry a large bulk-kinetic energy, 
tapping a large percentage of the energy available in the 
shock. A recent model of SQ, involving the collision be- 
tween two inhomogeneous gas flows, describes H2 forma- 
tion out of the multiphase, shocked gas, and an efficient 
cooling channel for high-speed shocks as an alternative 
to X-ray emission ([Guillard et al]l2009h 

Since the SQ detection, several other systems exhibit- 
ing similar l y strong H? emission have been discovered. 
lOgle et al.l (|2007l 120091 ) find that a large-subset of the 
local 3CR radio galaxies have extremely dominant MIR 
rotational H2 lines, often seen against a very weak ther- 
mal continuum. The low AGN and star formation power 
are insufficient to drive the MIR H2 emission. Mechan- 
ical heating driven by the radio jet interaction with the 



host galaxy ISM is the favoured mechanism. In addition, 
total H2 luminosities in the range 10 41 — 10 43 ergs s _1 
have bee n detected in some central cluster galaxies out to 
z ~ 0.3 ([Egami et al.ll2006l G. de Messieres - University 
of Virgin ia, Private commu nication) and in filaments in 
clusters ([Hatch et al.1 [2005) . The study of nearby pro- 
totypes may provide valuable insight into the nature of 
these more distant systems. The large scale (~ 30 kpc) 
of the SQ shock is well-suited for such a study. 

In th is Paper I, we extend t he single pointing observa- 
tions of lAppleton et all (|2006f ) to full spectral maps of SQ 
using the IRS instrument on the Spitzer Space Telescope, 
hereafter Spitzer. In Paper II, we will present detailed 
2-dimensional excitation maps of the H2 emission across 
the face of the X-ray emitting shock, and compare them 
with models. In addition, several other papers are being 
prepared by our team which will discuss the relationship 
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between the UV/X-ray emission and emission from dust. 

In Section 2, we present our observations and data re- 
duction methods. In Section 3 we discuss the mapping 
results for various detected lines. In Sections 4, 5 and 6 
we present our discussion on the properties of the system 
and in Section 7 the implication for galaxy formation. In 
Section 8 we present our conclusions. Additional mate- 
rial is included as appendices, with a discussion on NGC 
7319 in Appendix A, and a reanalysis of the high resolu- 
tion MIR spectrum of the shock, as well as renanalysis of 
the X-ray data presented in Appendix B and C, respec- 
tively. 

Throughout this paper, we assume a systemic velocity 
of v = 6600 km s -1 for the group, corresponding to a 
distance of 94 Mpc with H Q = 70kms~ 1 Mpc~ 1 . 

2. OBSERVATIONS AND DATA REDUCTION 

Mid-infrared spectroscopy of the shock region in SQ 
was obtained using the IRS instrument (jHouck et al.l 
2004) onboard Spitzer. Observations were done in low- 
resolution mapping mode, using the Short-Low (R ~ 
60 - 127; 5.2 - 14.5 fim) and Long-Low (R ~ 57 - 126; 
14.0 — 38.0 /Ltm) modules and taken on January 11 2008 
and December 10 2007, respectively. Figure 1 displays 
the outline of the areas observed superimposed on a com- 
posite image of the group. 

The SL spectral mapping consists of two separate, par- 
tially overlapping maps, centered North and South on 
the X-ray emission associated with the shock. The map 
was constructed with 23 steps of 2.8"(0.75 x slit width) 
perpendicular to the slit and one parallel step of 7.2". 
Observations consisted of 60 s integrations with 5 cycles 
per step. 

The LL module was used to map an area of ~ 2.8' x 3.2' 
using 21 steps of 8.0" (0.75 x slit width) perpendicular 
to the slit and a parallel step of 24.0". An integration 
time of 120 s was used with 3 cycles per step. 

Primary data reductions were done by the Spitzer Sci- 
ence Center (SSC) pipeline, version S17.0.4 and S17.2.0 
for SL and LL respectively, which performs standard re- 
ductions such as ramp fitting, dark current subtraction 
and flat-fielding. Background subtraction for LL data 
was performed by subtracting dedicated off-source ob- 
servations, with the same observing mode, taken shortly 
after the mapping sequence. In the case of SL, where for 
scheduling reasons the dedicated "off" observations were 
too far away in time to be optimal, backgrounds were 
generated from observations at the periphery of the map 
that contain no spectral line signatures. 

Examination of the pipeline products showed that the 
stray-light-corrected images (that account for potential 
spillover from the peak-up arrays onto the SL1/2 spec- 
tral apertures) contained wavelength-dependent over- 
corrections in some of the images, seemingly due to a 
high background of cosmic rays during the SL portion of 
the mapping. To negate this effect, the alternative flat- 
fielded images (also available from the science pipeline) 
that are uncorrected for crosstalk and straylight removal 
were used. It was determined that after background sub- 
traction, any uncorrected stray light in the SL spectral 
extraction areas was unmeasurable at the < 2 — 3% level. 
Thus stray-light correction was unneccesary, and the re- 
sulting Basic Calibrated Data (BCD) images used were 
of better quality than the standard BCDs. 



For all modules, individual frames for each pointing 
were median-combined and obvious "bad" pixels were 
replaced using customised software that allows for man- 
ual "average" replacement. The spectra were assembled 
into spectral cu bes for each modu le using the software 
tool, CUBISM ([Smith et al.ll2007l) . Further bad pixel 
removal was performed within CUBISM. Spectral maps 
were generated by making continuum maps on either side 
of a feature and subtracting the averaged continuum map 
from the relevant emission line map. One dimensional 
spectra were further extracted from the data cubes us- 
ing matched apertures. 

Broad-band images at 16 and 24/zm were obtained 
with the Spitzer IRS blue Peak-up Imager (PUI) and 
the MIPS instrument. The PUI was obtained in a 5 x 5 
map with 3 cycles of 30 s du ration each on 2007 December 
10. The MIPS instrument ([Rieke et al. Il2004h on Spitzer 
obtained 24/im imaging of SQ on 2008 July 29, achiev- 
ing a spatial resolution of ~6". Primary data reduction 
was done by the Spitzer Science Center (SSC) science 
pipeline (version S18.0.2) run through the MOPEX soft- 
ware, and for the MIPS image a smooth 2D polynomial 
background was further removed to correct for a large- 
scale background gradient. Spitzer IRAC 3.6, 4.5, 5.8 
and 8.0 /Ltm data of SQ (P.I. J.R. Houck) were obtained 
from the SSC archive; these were reduced using science 
pipeline version SI 8. 0.2. The final mosaics have a pixel 
scale of 0.61". 

3. RESULTS 

3.1. The Molecular Hydrogen Distribution 

The pure rotational transitions of molecular hydrogen 
can be excited by several mechanisms. These include 
FUV (Far Ultraviolet) induced pumping, and possible 
additional collisional heating, of the H2 in photodissoci- 
ation regions associated with star formation (e. g. Black 
& van Dishoeck 1987, Hollenbach & Tielens 1997), hard 
X-rays penetrating and heating regions within molecular 
clouds, which in turn excite H2 via collisions with elec- 
trons or hydrogen atoms (Lepp & McCray 1983; Drainc 
& Woods 1992) and finally collisional excitation of H2 
due to acceleration produced by shocks (e.g. Shull & 
Hollenbach 1978). The pure rotational MIR line ratios 
are not especially good diagnostics for distinguishing be- 
tween these mechanisms since all three mechanisms dis- 
cussed can lead to well thermalized level distributions of 
lower-level rotational states. The rotational H2 emission 
lines do, however, allow us to trace gas at different tem- 
peratures and compare with model predictions (this will 
be the main emphasis of Paper II). Higher level transi- 
tions 0-0 S(3)-S(5) tend to trace warmer gas, whereas the 
S(0) and S(l) lines are sensitive to the "coolest" warm 
H 2 . 

Although the line ratios themselves cannot be used 
directly as shock diagnostics, in Stephan's Quintet the 
distribution of large-scale X-ray and radio emission, plus 
optical emission line diagnostics, provide strong evidence 
that the giant filament seen in Figure 1 is the result of 
a strong shock. In Appleton et al. (2006) this fact was 
used to reveal the surprising association of detected H 2 
emission with the shock. However, in this paper we can 
make a more definitive association of the emission with 
the shock by means of spectral maps. 
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FlG. 2 . — Th e large-scale distribution of warm H2 in Stephan's Quintet with each detected 0-0 transition overlaid on an .R-band image from 
IXu et all d!999D . The grey box indicates the limit of the SL module coverage. Contour levels are as follows (all in units of MJy/sr): (a) 0.1, 0.14. 
0.19, 0.28, 0.32, 0.37, 0.41, 0.46, 0.5, (b) 0.3, 0.53, 0.75, 0.98, 1.2, 1.43, 1.65, 1.88, 2.1, (c) 0.3, 0.4, 0.5, 0.6, 0.7, (d) 0.25, 0.42, 0.58, 0.75, 0.92, 1.08, 
1.25, (e) 0.11, 0.19, 0.27, 0.36, 0.44, 0.52, 0.6 and (f) 0.3, 0.46, 0.61, 0.77. 0.93, 1.09, 1.24, 1.4. The SL coverage for (f) is larger due to incorporating 
the "parallel mode" SL2 slit observations from the "off-target" positions when SL1 mapping observations were the primary observing mode. 
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The spectral cubes were used to extract maps of all 
the pure rotational emission lines of molecular hydro- 
gen that were detected, namely the 0-0 S(0) 28.22/im, 
S(l) 17.03/im, S(2) 12.28/mi, S(3)9.66/im, S(4) 8.03/Ltm 
and S(5) 6.91/im lines. Specific intensity contour maps of 
these lines are presented i n Figure 2 overl aid on a i?-band 
image of the region from IXu et all (fl999h . The S(0) and 
S(l) lines were mapped by the LL modules, while S(2) - 
S(5) transitions were mapped by the SL modules of the 
IRS. As indicated in Figure 1, the SL observations were 
concentrated on the main shock to provide high signal to 
noise (S/N) measurements there. As a result, these maps 
do not fully cover SQ-A or NGC 7319. We note that the 
S(4) line at 8.03/xm (Fig. 2e) is faint, and also suffers 
from contamination from the PAH bands at 7.7/im and 
8.6/im. 

The contours indicate powerful, widespread emission 
running North-South along the shock ridge (see Fig. 1). 
In addition we detect strong emission from the star form- 
ing region, SQ-A, as well as associated with NGC 7319. 
We discuss further details of the AGN-like MIR emission 
lines from NGC 7319 in Appendix A. 

Figures 2a-d reveal a new H2 structure running east- 
ward from the "main" shock ridge. In what follows, we 
refer to this feature as the "brid ge". This structure is 
observed fa intly in the Chandra jTrinchieri et al.l 120031) 
and XMM (jTrinchieri et al.ll2005l) X-rav images and de- 
tected as faint Ha emission by IXu et al.l (|2003l ). but is 
not strong in radio continuum images. 

As is evident in Figure 2, there is distinct variation in 
the distribution of the warm H2 emission. The brightest 
0-0 S(0) emission (Fig. 2a) appears to be concentrated 
towards the North of the shock, whereas the 0-0 S(l) 
transition emission (Fig. 2b) appears more concentrated 
towards the center. The S(0) and S(l) maps demonstrate 
that the H2 in the bridge terminates in a large clump a 
few arcsecs west of the nucleus of the Seyfert 2 galaxy 
NGC 7319, and in a small detour to the North (especially 
in the S(l) map which has the highest S/N). 

The S(2) through S(5) lines clearly show that the warm 
H2 emission breaks into clumps in the shock. Despite the 
limited coverage compared to the LL mapping, the base 
of the "bridge" is visible and SQ-A is partially covered. 
SQ-A is fully covered by the SL2 module (because of 
fortuitous "off observation" coverage) and hence the S(5) 
emission line reveals that in SQ-A the H2 emission is also 
clumpy. 

The molecular-line maps provide considerable informa- 
tion about the excitation of the H2 along and across the 
X-ray shock, but these discussions will be deferred to 
a full 2-D modelling of the H2 excitation in Paper II 
(Appleton et al., in preparation). Instead we shall limit 
ourselves to global properties of the H2 here. In Section 
2|we shall present spectra of some selected regions of the 
emission and discuss a global excitation diagram for the 
shock. 

To further demonstrate the close connection between 
the H2 emission and the main global shock-wave in SQ, 
we now consider the distribution of warm H2 in relation 
to the X-ray and radio emission. Figure [3^, and b show 
the S(0) and S(3) contours overlaid on an XMM-Newton 
X-ray image (|Trinchieri et all [2005) of Stephan's Uum- 
tet. The warm molecular hydrogen is distributed along 
the length of the main North-South (NS) X-ray shock 



ridge and along the "bridge" , demonstrating the remark- 
able projected coexistence of of hot X-ray plasma (10 6 < 
T < 10 7 K) and warm H 2 (10 2 < T< 10 3 K). Although 
the H2 appears to follow the X-ray, there are subtle dif- 
ferences. The cooler S(0) line has emission concentrated 
to the North and follows the X-ray less closely compared 
to the warmer S(3) line. The S(3) line exhibits a clear 
correspondence to the X-ray, notably at the center of the 
shock, where we find peaks at both wavelengths. Thus 
the region of greatest shock heating, as traced by the 
stronger X-ray emission, appears to correspond to the 
higher- J H2 transitions, perhaps implying a causal con- 
nection. The intergalactic star formation region SQ-A, 
is essentially absent in X-ray emission, as seen in Figure 
[3^,, but is strongly detected in H2. 

We find a similar picture in the radio continuum (Fig. 
[3J: & d) with the S(0) line demonstrating correspon- 
dence with the main shock, but dominated by emission 
in the North where we observe less powerful radio emis- 
sion. The S(3) (and S(2)-see Fig. 2c) line presents 
a much tighter correlation with regions of the shock 
that are more radio luminous than the lower-J transi- 
tions. The radio emission is quite likely sensitive to the 
most compressed regions of the shock w here cosmic ray 
parti c les are accelerated m ore strongly (jHelou fe Bicavl 
Il993t I Appleton etHI 119991) . whereas the brightest X- 
ray pa tches are likely due to the fastest regions of the 
shock (jGuillard et al.l [2009). As already mentioned the 
"bridge" emission is detected only faintly in the X-ray 
and is weak or absent at radio continuum wavelengths. 
This noticeable difference compared to the main shock 
likely implies that the conditions that give rise to strong 
synchrotron emission in the main shock are absent in the 
bridge. 

3.2. Shock-related Fine- structure Line Emission 

Emission from fine-structure lines provide key di- 
agnostics that trace the interplay between the vari- 
ous constituents of the shocked interstellar medium. 
In lAppleton et al.l (|2006D the spectra were limited to 
the very core of the shock and only weak emission 
was detected from all but two metal lines, namely 
[Neii]12.81/zm and [Si ii]34.82/mi (this data has been 
reanalysed and is presented in Appendix B). In this 
section we discuss the spatial distribution of emission 
from the [Feii]25.99/im, [O iv]25.89/im, [S m]33.48/im, 
[Siii]34.82/zm, [Neil] 12. 81 /urn and [Nelll]15.56 //m lines. 

In Figure [5] we present the specific intensity con- 
tours of the [Feii]25.99^m and [O iv]25.89/im blend, 
[S ni]33.48/zm and [Si n]34.82^m emission lines. 

Given the low spectral resolution of the SL and LL 
modules of Spitzer, we cannot distinguish between emis- 
sion from [Fe n]25.99/im and [O iv]25.89^tm. However, 
except in the direction of the Seyfert II galaxy NGC 7319, 
the emission near 26/im seen in Figure [3Ji is likely to 
be pure [Fe 11] with little contamination from [O iv] as 
there is no evidence from the spectra of high-excitation 
emission from the intragroup medium in SQ. For exam- 
ple, [Oiv] has an excitation pot ential of 56 eV (com- 
pared to 7.9 eV of [Fe 11]), and vet lAppleton et all (|2OO60 
have shown that the [Nein]/[Neii] ratio is low suggesting 
low-excitation conditions for the ions in the shock, fur- 
ther supported in Section [5[ Assuming that, apart from 
towards NGC 7319, [Fen] dominates the [Feii] + [Oiv] 
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Image: 20 cm 



Fig. 3.— X-ray image from XMM-Newton ijTrinchieri et al.ll2005l ) in the 1.5-2.5 keV band with (a) H 2 S(0 ) contours and (b) 
H2 S(3) contours overlaid. Also shown are 1.4 GHz radio continuum images from the VLA l|Xu et al.1 120031 ) with (c) H 2 S(0) 
contours and (d) H2 S(3) contours overlaid. Contours levels are as in Fig [2] 



complex, we detect faint emission from [Fe n]25.99/im at 
the location of the center of the shock (as defined by the 
X-ray "hotspot" in Fig. [3^). The energetic requirements 
for shocks to produce strong [Fen] emission are usually 
present in J shocks whi le the ion abundance in C s hocks 
are low in comparison (|Hollenbach fe McKed fl989K We 
will discuss the production of [Fe 11] in the shock in Sec- 
tion 6.2. 

The [S in]33.48/im distribution is presented in Figure 
This f ine structure lin e acts as a strong tracer of H 11 
regions (jDale e t al. 200f|) and we observe emission from 
SQ-A and from other regions of star formation in the 
South (see Section 13. 3|) . In the primary shock region the 
distributions of [Fen] and [S ill] are anti-correlated. 

In strong contrast to the weak [Fe ulemission, there is 
copious [Si ii]34.82/im emission (Fig. which follows 
the S(l) distribution closely with respect to the primary 
shock, as mapped by X-ray emission (Fig. [Hi). 

Although [Si 11] is commonly found in normal H II re- 
gions, we will demonstrate below that, apart from in SQ- 
A, the strong silicon emission does not correlate with re- 
gions of strong PAH emission (tracing star formation) in 
SQ, but instead closely follows the H2 and X-ray (shock) 
distributions. [Sin] acts as an efficient coolant of X-ray- 



irradiated gas and is predicted to be one of the top four 
coolin g lines under these circumstances (Maloney et al. 
I1996T) . We will discuss the excitation of [Si 11] in the shock 
in Section 6.2. 

[Neii]12.81/im (with an ionisation potential of 21.6 eV) 
is also represented in the shock, as shown in Figure 
although it is also emitted from some Hn regions asso- 
ciated with star formation, such as SQ-A and the star- 
forming region South of NGC 7318b (see Section T3.3p . 
The [Nein]15.56^m contours (with an ionisation poten- 
tial of 41 eV) are shown in Figure [SJd and are associated 
with excitation in the shock, as well as from regions of 
star formation (see Section 13 . 3|) . The higher ionisation 
line of [Nein] is much weaker in the shock ridge, but 
regions of emission correspond closely to peaks seen for 
the H2 S(3) line. The [Neil] emission, however, exhibits 
clearly extended emission with regions of greatest lumi- 
nosity matching those seen for the H2 S(3) line. The 
[Ne 11] emission suggests excitation from the shock, with 
[Ne in] found at the location of the center of the shock, 
similar to what is observed for [Fe 11] . 

3.3. Star Formation Tracers and Dust Emission in 
Stephan's Quintet 
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Fig. 4. — Fine structure emission ([Fe n]25.99/im and [O iv]25.89/im blend, [S ni]33.48/im and [Si ll]34.82/im) compared to the 
distribution of Eh S(l) emission. Also shown are [Si ll]34.82/im contours compared to X-ray 0.5-1.5 keV emission. Contour 
levels are as follows (all in units of MJy/sr): (a) 0.08, 0.1, 0.13, 0.15, 0.165, 0.18, 0.2, 0.22, 0.23, 0.25, (b) 0.2, 0.33, 0.47, 0.6, 
0.73, 0.87, 1.0 with (c) and (d) having levels of 0.5, 0.66, 0.83, 0.99, 1.15, 1.31, 1.48, 1.6, 1.8. These levels are chosen to best 
show intensity of emission associated with the shock and do not reflect the full intensity range for NGC 7319. 



Previous obs ervations and spectroscopy of SQ 
(|Xu et al.l 120051 ) have determined that there are some 
regions of star formation associated with the spiral arms 
of the intruder galaxy NGC 7318b. We discuss in this 
section how these regions, which have a different spatial 
distribution from the shocked gas, are correlated with 
the PAH emission we detect in the IRS spectra. 

In Figure [6^ there is a strong correlation between the 
11.3/tm PAH distribution (from the IR S cube) superim - 
posed on a Near-UV image from Galex (|Xu et al.ll2005h . 
which maps the UV emission from hot stars associated 
with weak star formation from the system. This correla- 
tion suggests that the PAH molecules are excited by star 
formation. A similar close correlation is shown in Figure 
\Ejp> where we overlay the 11.3/xm contours on the IRAC 
8/im band, which is dominated by the 7.7 and 8.6 //m fea- 
ture. It is noticeable, however, that regions with strong 
11.3y«m emission in the shock, do not appear similarly 
strong at 8/im. This point will be addressed in Section 



6.4. 

We compare the 11.3/xm PAH map to the distribution 
of warm dust in Figure Et, using the MIPS 24/im map 
of SQ. Again, it is clear that there is a good correlation 
between the PAH emission and the thermal dust, most of 
which seems only poorly correlated with the shock ridge. 
The lack of cons picuous st ar formation in the ridge was 
also observed bv lXu et all (|2003D . 

The main point we emphasize here is that the dust, 
PAH and UV emission appears to be associated with pre- 
viously known star formation regions and no additional 
star formation is observed in the shock; this can be seen 
in Figure [5Ji showing the H2 0-0 S(l) emission overlaid 
on an IRAC 8/im image. There is little correspondence 
between the H2 emission in the shock and 8//m (hot dust 
plus PAH) image. This is important because it implies 
that there is very little triggered star formation in the 
molecular gas associated with the shock-excited H2. 

Figure demonstrates that there is only faint dust 
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(a) [Ne II] contours 



(b) [Nc 



contour 




Image: H 0-0 S(1; 



Fig. 5. — [Nen]12.81^m and [Nenl]15.56/^m emission compared to the distribution of H2 S(l) emission. Contour levels are as 
follows (in units of MJy/sr): (a) 0.22, 0.4, 0.57, 0.75, 0.92, 1.1 and (b) 0.06, 0.08, 0.1, 0.13, 0.16, 0.18, 0.2, 0.23, 0.25. 



emission at 24/xm from the shock ridge. The presence 
of dust in the shock is required in the model of Guillard 
et al. (2009) to explain the formation of H2 behind the 
shock and we do observe evidence of depletion onto dust 
grains (see Section 6.2). Thus the faint 24/xm emission 
could be the result of destruction of Very Small Grains 
(VSGs), with only larger grains surviving, or indicate 
that the grains are cold and radiating more strongly at 
longer wavelengths, where Spitzer has the least spatial 
resolution. A more detailed description of the dust and 
faint PAH emission in the SQ group (including results 
from MIPS 70/im imaging) is discussed in separate pa- 
pers (Guillard et al. 2010.; Natale et al. 2009, in prepa- 
ration). Guillard et al. (2010) show that the IR emission 
in the shock is faint due to dust being heated by a rel- 
atively low intensity UV radiation field and determine a 
Galactic PAH/VSG abundance ratio in this region. 

A more complete understanding of the likely existence 
of cool dust in the shock will require higher angular res- 
olution and a broader wavelength coverage than that 
achieved by Spitzer. 

Despite the faintness of emission from the main shock, 
the 24/im map presents a new result, which was less ob- 
vious in previous studies, namely that the dominant re- 
gions of star formation in SQ lie not in the galaxies them- 
selves, but in two strikingly powerful, almost symmetri- 
cally disposed regions at either end of the shock. The 
region to the North is the well studied SQ-A, but the 
region to the South (which we refer to as 7318b-south) 
is also very powerful and both regions lie at the ends of 
the shock, as defined by the H2 distribution. This may 
not be a coincidence, and we will discuss this further in 
Section 6.1. 

The [Neil]12.81/mi, [Nem]15.56/mi and [S ni]33.48/xm 
fine-structure lines, as mentioned above, are also tracers 
of star formation in SQ as these lines are often associated 
with H 11 regions. The [Ne 11] 12. 81/im emission appears 
to follow both the H2 and the star formation regions (see 
Fig. [5^i), appearing more extended in the South than 
the corresponding H2 emission and flaring out where star 



formation regions, especially 7318b-south, are observed 
optically, and through PAH emission (see Fig. [S]). 

4. GLOBAL PROPERTIES OF THE MOLECULAR 
HYDROGEN IN SQ 

The warm molecular hydrogen in Stephan's Quintet 
follows the X-ray distribution in the main shock and 
in the "bridge" structure. This might suggest that the 
molecular hydrogen is excited directly by the X-ray heat- 
ing. However, we will show that the H2 emission exceeds 
by at least a factor of 3 the X-ray luminosity from the 
various shocked filaments, thus ruling out direct X-ray 
excitation from the shock. 

To measure the strength of the H2 emission, we extract 
spectra from various rectangular sub-areas of SQ which 
are defined in Figure [7] The spectra were extracted from 
CUBISM cubes built from each IRS module, and joined 
to make a continous spectrum-no scaling was necessary 
to join the spectra. 

Figures [7] a, b and c indicate the spectral extraction 
regions of the main NS shock, a sub-region of the shock 
and a characteristic part of the "bridge" . The shock 
sub-region is chosen to be just North of the center of 
the shock, avoiding regions contaminated by star forma- 
tion in the intruder. Spectra for these extractions are 
shown in Figure [H] a, b and c. All three spectra, share 
the common property that they are dominated by molec- 
ular hydrogen emission. The shock sub- region (Fig. [SJd), 
unlike Figure^, is less contaminated by the star form- 
ing regions discussed in the previous section. The MIR 
continuum of the main shock appears stronger than the 
shock sub-region indicating stronger emission from warm 
dust; this is likely the result of contamination from star- 
forming regions in the main shock extraction. The bridge 
exhibits a similarly weak continuum emission compared 
to the shock sub-region. 

Figure [BJa includes photometry from the IRAC bands, 
16/im Peak-up Image (PUI) and MIPS-24/im image su- 
perimposed on the IRS spectrum. These are useful to 
probe conditions in the shock, in particular star forma- 
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(a) 11.3 urn PAH contours 




Fig. 6. — Dust as tracers of star formation. The 11.3/im PAH 
a 24/im image. By comparison the H2 0-0 S(l) contours show 
levels are as follows (in units of MJy/sr): (a), (b) and (c) are ( 



(b) 1 1 .3 um PAH contours 




contours overlaid on (a) the Near-UV, (b) an 8/im image and (c) 
little correlation with the 8fim image (as shown in d). Contour 
3.2, 0.35, 0.5, 0.65, 0.8. Contour levels for (d) are as in Fig. [2b. 




(a) 



(b) 




(c) 



(d) 



Fig. 7. — The rectangular sub-regions extracted from the main components of the SQ shock system overlaid on the H2 0-0 
S(l) line emission, a) The Main Shock region centered on 22 h 35 m 59.6 s , +33°58'05.7". b) A region in the main shock centered 
on 22 h 35 m 59.97 s , +33°58'23.3"and chosen to have limited contamination from star forming regions in SQ. c) Extraction of the 
H 2 "bridge" feature centered on 22 h 36 m 0.14 s , +33°58'23.3". d) An extraction of the SQ-A star forming region centered on 
22 h 35 m 58.85 s , +33°58'50.4". 
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FlG. 8. — Extracted Spectra (low resolution) of the regions shown in Figurc[7] combining all IRS modules with coverage. The key emission features 
arc labeled with measurements and sizes of extraction areas given in Table [2] which provides the H2 line fluxes, and Table [3] listing the metal line 
measurements. The shock sub-region (shown in b) includes the IRAC, 16/j-m PUI and 24/xm MIPS measurements for the same area. The 16/xm flux 
is 0.939±0.047 mjy, which includes emission from the 17.03/xm H2 line covered by the filter, while the 24/im flux is 0.40S±0.046 mjy; photometry 
was performed using a circular aperture with no aperture correction applied to compensate for contamination from the surrounding emission. The 
corresponding IRAC fluxes are 1.362±0.068, 0.885±0.045, 0.659±0.036 and 1.434±0.073 mjy for the 3.6, 4.5, 5.8 and 8.0 /xm bands, respectively. 
The extraction of SQ-A (shown in d) also has coverage from the SL2 module due to the orientation of the IRS during the SL1 observation. 
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tion (see Section 6.4). Stellar light in the extraction area, 
from an extended spiral arm of NGC 7318b, produces 
contamination of the shock spectrum, visible as contin- 
uum emission shortwards of 6/im, in both Figure 8a and 
b. 

There exists a striking similarity between the mid-IR 
spectrum of all three regions, showing powerful H2 lines 
and low excitation weak emission from fine structure 
lines. Also the PAH emission observed in the spectra of 
the shock regions, and in the region of the "bridge" in the 
IRAC 8/im image (Fig. [6b), appears wea k. This confirms 
that t hese properties, observed in the Applcto n et alJ 
(2006) observations of the shock core, extend to both 
the full extent of the main shock and the "bridge" . This, 
and the fact that the "bridge" has similar X-ray prop- 
erties to the main shock (see later) suggests that the 
bridge is a "scaled-down" version of the main shock. The 
weaker radio continuum emission at this location is sig- 
nificant. One possibility, that the "bridge" is older than 
the main shock, and the cosmic rays compressed in it 
have diffused away, will be discussed further in Paper 
II. Thus the new IRS observations seem to suggest that 
more than one large-scale group-wide shock is present in 
the group. This could be the result of previous tidal in- 
teractions and imply multiple shock heating events have 
taken place in SQ, consisten t with what is seen in the 
X-ray (|Q'Sullivan et~alll2559h . 

Figures [31 and [8ji present the extraction region and 
spectrum of SQ-A, the extragalactic star-forming region. 
In this case, although H2 lines are still strong, a rising 
continuum and an increase in the strength of the metal 
lines relative to the H2 is consistent with a spectrum that 
is increasingly dominated by star formation - a result 
which is already kno wn from previous optical observa- 
tions (|Xu et alfeool . 

Line fluxes for all the H2 and metal lines in the spectra 
discussed above are presented in Tableland Table [3J re- 
spectively. We estimate the luminosity emitted from the 
H2 lines in the main aperture shown in Figure [8^,. The 
emission from the 0-0 S(l) line alone can be calculated 
from Tabled] (for D = 94 Mpc) to be 2.3 xlO 41 ergs" 1 . 
Summing the emission measured in the observed lines 
for the main shock (0-0 S(0) through S(5) lines) and in- 
cluding an extra 28% emission from unobserved lines (see 
model fit to excitation diagram below) , yields a total H2 
line luminosity from the main shock of 9.7 x 10 41 ergs -1 . 
This phenomenal power in the molecular hydrogen lines 
dwarfs by a factor of ten the next brightest mid-IR line, 
which is [Si n]34.82/im with a line luminosity of Lgin = 
0.85 xlO 41 ergs" 1 . 

Figure [H] presents the excitation diagram of the low- J 
0-0 H2 transitions for the main shock extraction. The 
points are well fit by a model including three temper- 
ature components (Ti = 158±6K, T 2 = 412±4K, and 
T 3 = 1500±200K). It is likely that in reality, many dif- 
ferent temperature components are present in the shock, 
and the three-temperature fit is only an approximation. 
However, it does allow us to provide an estimate of the 
total mass of warm H2 of 5.0±lxl0 8 M Q . Temperature 
T2 is more uncertain than formally represented by the fit 
because it depends on the value of the S(4) flux, which 
may be systematically too low due to PAH contamination 
(see Fig. [§]). In Paper II we will present a more com- 



plete two-dimensional map of the excitation of the H2 in 
SQ and explore variations in the shape of the excitation 
diagram along and across the shock in more detail. 

Our observations have shown that H2 is the dominant 
line coolant in the MIR. However, how does it compare 
with the most important co olant in high-sp e ed sh ocks- 
namely the X-ray emission? lAppleton et al.l (|2006l ) sug- 
gested that the H2 emission was stronger than the X-ray 
emission at the shock center. This can now be evaluated 
over much of the inner SQ group. 

We present a complete reanalysis of the XMM-Newton 
observations of SQ using the latest calibrations (see 
Appendix C for full details) in order to determine 
the fluxes and luminosities of the X-ray emission to 
match our spectral extractions. The results indicate 
the striking dominance of the H2 line luminosities com- 
pared with the X-ray emission from the same regions. 
For the main shock, the X-ray "Bolometric" flux of 
Lx(o.ooi-iokcV) = 2-8 x 10 _13 erg s _1 cm -2 corresponds 
to 2.95 x 10 41 erg s" 1 , or L(H 2 )/L x(0 . 001 _ 10keV) = 2.9. 
This is a lower limit since we have not attempted to 
remove the contribution to the main shock aperture of 
an extended group-wide X-ray component upon which 
the emission from the shock lies. Therefore it is likely 
that the H2 line luminosity dominates over the main 
shock X-ray gas by a factor > 3. Similar calculations 
can be done for the other regions for which H2 spec- 
tra have been extracted. For example, in the "bridge" 
region, which we have already indicated has many of 
the same characteristics as the main shock, we find 
L(H2)/Lx(o.ooi-iokeV) = 2.5. These values demonstrate 
that throughout the extended regions of SQ, the molecu- 
lar hydrogen cooling pathway dominates over the X-ray 
in this shocked system. This is a very significant result, 
upturning the traditional view that X-ray emission al- 
ways dominates cooling in the later stages of evolution 
in compact groups of galaxies. 

5. EMISSION LINE DIAGNOSTICS 

The fine-structure flux ratios (see Tabled can be used 
to probe the conditions within the extracted regions of 
SQ. The [Siii]34.82/im/[Sm]33.48/xm ratio provides an 
indication of the sources of excitation within the sys- 
tem. As mentioned previously, [S in] is mainly a tracer 
of Hn regions, whereas enhanced [Sin] emission can be 
generated via several mechanisms, including thermal ex- 
citation by X-rays (XDRs), or in shocks. In the main 
shock, the [Si n] 34.82^m/[S IIl]33.48/zm ratio of - 4.59 is 
high compared to, for example, both normal galaxies (~ 
1.2) a nd AGN (~ 2.9) in the SINGS sample (jfaale et all 
2006). However, this large aperture is contaminated 
by star formation emission from SQ-A and the intruder 
galaxy. A better measure is given by the smaller shock 
sub-region, where the [Si n] 34.82^m/[S iii]33.48/xm ra- 
tio is ~ 9.5. Thus it is clear that the [Sin] emission 
is well outside the normal range of values, even for lo- 
cal well-studied AGN. Using the upper limit found for 
[S Hi] 33.48/xm in the "bridge" structure, we find a ra- 
tio of > 5 - again values well outside the range of nor- 
mal galaxy disk emission. Indeed, these high values 
are typical of galactic supernova re mnants where shock 
excitation is well determined (e.g. iNeufeld et al.l 120071 ; 
IHewitt et all 120091 . We will argue in the next section 
that silicon is being ionised in regions experiencing fast 
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Fig. 9. — Excitation Diagram of the Main Shock region in Stephan's Quintet. The ordinate (N/g) represents the column density 
divided by the statistical weight of the rotational state. 



shocks 100 < V s < 300 km s" 1 and depleted onto dust 
grains. The [Fen] 26.0/mi/[Siii] 34.82^m ratio for the 
main shock and sub-region A 0.12) i s in agreement with 
values found by iNeufeld et a l. ( 2007) for their sample of 
SNR. 

For the extragalactic star forming region SQ-A, 
[Sin] 34.82//m/[S m] 33.48/mi is ~1.49, only slightly 
higher than the average of ~ 1.2 found for sta r- forming 
regions m the SINGS sample pale et al.l l2006): another 
indication that SQ-A is dominated by star formation. 

In star-forming galaxies the 

[Nelll] 15.56^m/[Neii] 12.81/xm ratio can be used 
as a measure of the hardness of the radiation field as it 
is sensitive to the effective temperature of the ionising 
sources. In the main shock, we find a value of ~ 0.36 
which would be considered typical compared to those 
found in starburst s y stems, which range from ~ 0.03 — 15 
(jVerma et al.1 120031: iBernard-Salas et all l2009h . and in 

supernovae remnants ranging fro m ~ 0.07 — 1.24 

(|Neufeld et alj [20071 : iHewitt et all [20091) . The shock 
subregion has a ratio of only ~ 0.14 indicating a lower 
intensity radiation field North of the shock center. 
However, it is clear from its spatial distribution relative 
to the 8/mi and H2 emission, that most of the neon is 
not originating from star formation, and so shocks are 
an obvious source of excitation. 

The [Neil] 12.81/xm/[Nem] 15.56^m ratio can be used 



to estimate the shock velocity us i ng th e MAPPINGS 
shock model library of lAllen et al.l (|2008l ). In the main 
shock this ratio corresponds to shock speeds of be- 
tween 100 and 300 km s -1 (using preshock densities of 
0.01 < n < 100 cm -3 and magnetic parameter B/n 5 
= 1 and 3.23 - the nominal equipartition value 11 ). We 
do, however, have contamination from starforming re- 
gions in the main shock and cannot disentangle this 
emission from that produced in the shock. To ad- 
dress this we use the [Neil] 12.81/xm and [Nein] 15.56/im 
maps to mask areas associated with star formation 
(see Section 3.3) and determine a lower limit for the 
[Nell] 12.81/xm/[Neiii] 15.56/xm ratio in the shock of 
~ 4.54. This value corresponds to a shock velocity of 
~ 150kmsA 

The average electron density is determined from the 
[Sm]18.71/xm/[Siii]33.48^m (two lines of the same ion- 
isation state) ratio. For SQ-A we find a ratio of 
~0.56, and in the main shock ~0.41. This cor- 
responds to an electron dens ity of 100 — 200 cm -3 
(|Martm -Hernandez et al. 2002T) for both, i.e. in the low- 
density limit for this diagnostic (|Smith et al.l l2009). 

6. THE GROUP- WIDE SHOCKS 

6.1. Origin of the H2 and X-ray emission 

11 See Allen et al.(2008) 
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Our observations have shown that the molecular hy- 
drogen and X-ray emitting plasma appear to follow a 
similar distribution, and we have ruled out the possibil- 
ity that this is a consequence of X-rays heating the H2 , 
since the H 2 has the dominant luminosity. How then 
can we explain the similar distributions? Are these re- 
sults consistent with the hypothesis that the shock is 
formed where the intruding galaxy NGC 7318b collides 
with a pre-existing tidal filament of H 1 drawn out of 
NGC 731 9 in a previous interaction with another group 
member ([Moles et al.lll997t iTrmchieri et al.ll2003ft ? 

This basic mechanical picture appears plausible as can 
be seen in Figure[TT)lwhich shows that the H2 distribution 
"fills i n" the gap in the H i tidal tail as observed by the 
VLA ([Williams et al.ll2002ft . The implication is that the 
H 1 has been converted into both a hot X-ray component 
and a warm H2 component by the collision of the intruder 
with the now missing H I. 

Part of the puzzle of how this high-speed (V s ~ 700 — 
1000 kms -1 ) shock can lead to both X-ray and very 
str ong molecular li ne em ission is presented in a model 
by IGuillard et~aH (|2009ft . The high-speed collision of 
NGC 7318b with the Hi filament (assumed to be com- 
posed of a multiphase medium) leads to multiple shocks 
passing through and compressing denser clumps (which 
become dusty nucleation sites for H2 formation) as op- 
posed to the lower-density gas, which is shock- heated to 
X-ray tempertures. The H2 therefore forms in denser 
clouds experiencing slower shocks. Thus the coexistence 
of both hot X-ray gas, and cooler molecular material is 
a natural consequence of the multiphase medium of the 
pre-shocked material. 

Modeling of the H 2 excitation bv IGuillard et all (|2009l ) 
demonstrates that the emission can be reproduced by low 
velocity (~ 5 — 20kms~ 1 ) magnetohydrodynamic shocks 
within the dense (tih > 10 3 cm~ 3 ) H2 gas. The denser 
clouds survive long enough to be heated by turbulence 
in the hot-gas component, tapping into the large avail- 
able kinetic energy of the shock. This picture is con- 
sistent with both the broad H2 linewidth (870 km s -1 ) 
measure d in IRS high-r e soluti on spectrometer observa- 
tions of I Appleton et alJ (|2006ft . and the velocity center 
of the warm H2 (based on new IRS spectral calibrations 
- see Appendix B) which places the gas at intermedi- 
ate velocities between the intruder and the group IGM. 
Both these measurements are consistent with H2 being 
accelerated in a turbulent post-shocked layer. 

Intermediate pre-shock densities and post-shock tem- 
peratures result in regions of Hi and H11 that have 
cooled, but where th e dust content has been destroyed 
(|Guillard et alJ 12009ft . Pre-existing Giant Molecular 
Clouds (GMCs) which may have been embedded in the 
Hi gas, would be rapidly compressed and collapse 
quickly, t hus forming stars . This mechanism, proposed for 
SQ-A bv IXu et all (|2003ft . might also apply to 7318b- 
south. However, if this was the case, it would have to 
explain why two such GMCs happened to be positioned 
at the extreme ends of the current shock - an unlikely 
coincidence. More probable, however, is that the geom- 
etry of the shock somehow favors the collapse of clouds 
at the ends of the shock - perhaps in regions where the 
turbulent heating is less efficient. 



6.2. Origin of [Fell] and [Sill] emission 

As outlined above, the combination of emission de- 
tected in the shock region of SQ can be understood in 
terms of a spectrum of shock velocities. The fastest shock 
velocities (V s ~ 600 kms -1 ) are associated with the low- 
est density pre-shock regions and the post-shock X-ray 
emitting plasma. These are fast J shocks and represent 
a discontinous cha nge of hydrodynamic variab les and are 
often dissociative (|Hollenbach fe McKed fl989). C shocks 
have a broad transition region such that the transition 
from pre-shock to pos t-shock is continuou s and are usu- 
ally non-dissociative (|Draine et~aT]H983l ). The lowest 
velocity shocks associated with th e turbulent Ho emis - 
sion are ~ 5 — 2 kms -1 C shocks (jGuillard et al.ll2009ft . 
IXu et all ([2003ft find optical emission line ratios consis- 
tent with shock models that do not i nclude a radiative 
precursor (jDopita fc Sutherlandlfl995ft . 

The [Fe ii]25.99/im emission associated with the shock 
region in SQ is relatively weak, but coincides with the 
most energetic part of the shock as traced by the X- 
rays (Fig. [4|). We also detect abundant [Si ii]34.82/im 
emission associated with the main shock. 

Silicon and iron have very similar first and second ion- 
ization potentials. Their first ionization potentials (7.9 
and 8.15 eV for Fe and Si, respectively) are lower than 
that of hydrogen but their second ionization potential 
is higher (16.19 and 16.35 eV). The MIR [Fen] and 
[Sin] line emission observed from the SQ shock could 
thus arise from predominantly neutral, as well as ion- 
ized gas. We discuss the contribution from the ion- 
ized gas using the [Fen](25.99/mi)/[Nen](12.81/im) and 
[Si 11] (34.82/im)/[Ne 11] (12.81/im) line ratios. The high 
[Neil] (12.81/im)/ [Ne ill] (15. 56/mi) mid-IR line ratio (see 
Section 5) implies that [Neil] is the dominant ionization 
state of Ne in the SQ shock . Unlike Fe and Si, Ne is 
not much depleted on dust ([Simpson et al. 1998). Be- 
cause the [Neil] 12.81/im line has a high critical density 
(n e = 4.3 x 10 5 cm -3 , Ho and Keto 2007), the neon line 
strength scales with the emission measure of the ionised 
gas. 

The optical line emission from the SQ shock is dis- 
cussed in detail in Xu et al. (2003). The high val- 
ues of the [Oi](6300A) and [Nn](6584A) lines to Ha 
line ratios are evidence of shock ionization. The optical 
[Sn](6716A/673lA) line ratio as well as the MIR [S in] 
line ratio correspond to the low density limit (see Section 
5) and comparison with the shock models of Allen et al. 
(2008) constrain the pre-shock gas density to be about 
1 cm -3 or smaller. 

In Figure [11] we indicate the region (in grey) corre- 
sponding to the observed MIR ratios (using the upper 
and lower limits determined for the [Ne 11] 12.81/im line 
emission as discussed in Section 5). For comparison the 
expect ed emission from the shock models of I Allen et alJ 
( 2008) are displayed for a pre-shock gas density of 1 cm -3 
and two values of the magnetic parameter, B/n 5 = 0.5 
and 3. For clarity of the figure, only shock velocities 
from 100 to 300 kms -1 are used and higher velocities 
shocks which do not match the observed [Ne 11] /[Ne in] 
line emission ratio (see Section 5) are discarded. 

The observed [Fen](25.99/im)/[Nen](12.81/mi) and 
[Si 11] (34.82/xm)/ [Neil] (12.81/mi) ratios are both smaller 
than the shock values. The iron and silicon lines are 
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Fig, 10-— Sch ematic Diagram of the Hi distribution in Stephan's Quintet, from William s et al] (j2002D . and the H a emission 
ijXu et al.ll 19991 ). with the H 2 0-0 S(l) emission overlaid. The grey dots represent the locations of the core galaxies in the group 
as shown in Fig. [1] 



not dominant cooling lines of ionizing shocks. The gas 
abundances of these elements do not impact the thermal 
structure of the shock and the line intensities roughly 
scale with the gas phase abundances. We thus interpret 
the offset between the IRS observation and model values 
as evidence for Fe and Si depletion. We consider the mag- 
nitude of the depletions indicated by the arrow (~50% 
and ^60% for Fe and Si, respectively) in Figure QT] as 
lower limits since there could be a contribution to the 
[Fell] and [Sin] line emission from non-ionizing J -shocks 
into molecular gas (|Hollenbach fe McKed I1989T ) . Such 
shocks could also contribute to the [O i] (6300A) line emis- 
sion as discussed in Guillard et al. (2009). Forthcoming 
observations of the H2 ro- vibrational line emission in the 
near-IR should allow us to estimate whether they may be 
significant. A contribution from non-ionizing shocks to 
the [Fe 11] and [Si 11] line emission will raise the depletion 
of both Fe and Si as well as the Fe/Si depletion ratio 
because non-ionizing shocks do not produce [Neil] line 
emission and the [Feil](25.99/im)/[Sin](34.82/zm) line 
emission ratio in J-shocks within dense gas is larger than 
for shocks plotted in Figure QT] 

6.3. Comparison with Cold H% Distribution 

A key aspect to understanding the emission we observe 
in SQ, as well as a test of the proposed model of H2 ex- 
citation, is the amount and distribution of cold H2, as 
measured by CO. A reservoir of cold H2 associated with 
the warm gas observed in SQ would provide key insight 
into conditions within the intragroup medium. In Fig- 



ure [12] we sho w the BIMA CO (1-0) integrated intensity 
contours from lGao fe Xul (|2000f h shown with the H 2 S(0) 
and S(l) emission contours, overlaid on an optical image 
of SQ. These interferometric observations use the large 
primary beam of 110", centered on the shock, to deter- 
mine areas of high column densities in the group. The 
CO traces areas of known star formation seen at 24/im 
and in the NUV (see Section [3. 3D . not ably S Q-A and 
7318b-south. As observed bv lYun et all (|1997h . the CO 
around NGC 7319 is concentrated in two regions. The 
dominant complex is North of the nucleus residing in a 
dusty tidal feature. The nuclear CO is elongated per- 
pendicular to the stellar disk suggesting a deficiency of 
ongoing star formation in the disk. 

The CO distribution does not correlate with the lo- 
cation of warm H2 emission, particularly around NGC 
7319. The cold H2 complexes are clearly offset from the 
concentrations of warm H2. Even the S(0) emission line 
(see Fig. fl"2"k). which follows the coldest warm H 2 , does 
not have peaks corresponding to the strongest CO detec- 
tions. 

In a forthcoming paper (Guillard et al. 2010, in prepa- 
ration), we will report the recent detection of 12 CO(l- 
0) and (2-1) emission, associated with the warm H2 in 
the SQ shock, using the single-dish 30m IRAM telescope. 
These observations suggest that most of the CO emission 
in the shock has been missed by interferometers (because 
of the broad linewidth) and show that the CO emission 
is both present in the shock, and extends along the H2- 
emitting bridge and towards NGC 7319. 
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Fig. 11. — The line ratios measured in the Main SQ shock of the [Nell](12.8/im), [Fell](26.0/im), and [Sill](34.8/im) lines are 
compared with shock calculations from Allen et al. (2008). The model values cover a shock velocity range of 100 — 300 kms - 
with a gas density of 1cm -3 and two values of the magnetic parameter B/n°'° — 0.5 (green square) and 3 (red circles). The 
black symbols represent the 100 kms -1 values and the arrows indicate the direction of increasing velocity (in steps of 25 kms -1 ). 
The grey area indicates the range of line ratios, [Sill]/[Nell] and [Fell] /[Nell], determined from observations (as discussed in 
Section 5). The ratios are lower than the model values which we interpret as evidence for Si and Fe depletion onto dust grains. 




Fig. 12. — CO (1-0) contours (green ) from iGao fe Xul (gOOO) overliad on an R-band image l|Xu et al.lH999l ) of the group. 
Contour levels correspond to 2.4a, 3.0a, 4.0a, 5.0cr, 6.0a, 7.0a and 8<r (I.Oct ~0.025 Jy). We show the H2 contours (red) for the 
S(0) line (a) and the S(l) line (b) with contour levels as in Fig. [2^ and b respectively. 

The kinematics of the CO gas lying outside of star- highly disturbed with a broad linewidth in agreement 
forming regions, in the new observations, appears to be with the Appleton et al. (2006) interpretation that the 
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MIR H2 lines were intrinsically very broad and resolved 
by the high resolution module of IRS. The CO data also 
agrees with our re-analysis of the Appleton et al. (2006) 
data (see Appendix B) using more recent and reliable 
wavelength calibration, which places the bulk of the H 2 
gas at velocities intermediate to that of the intruder and 
the group - supporting the idea that the H 2 gas is accel- 
erated in the shock. 

6.4. Turbulent Suppression of Star Formation in the H2 

filaments 

Section l3~3l discussed the star forming regions observed 
in SQ and noted that there was very little evidence for 
star formation in the shock associated with the warm 
H2 emission. The main shock contains a large quan- 
tity of warm molecular hydrogen (5.0 xl0 8 M Q ) provid- 
ing a reservoir of fuel for star formation once it cools 
(see Section [4j. We investigate star formation in the 
warm ^-dominated medium by considering the shock 
sub-region, chosen to avoid star forming regions in the 
intruder galaxy, but likely still subject to some contami- 
nation from these regions. 

PAH emission is a classical tracer of star formation, 
but the molecules are fragile and easily destroyed in hard 
radiation fields. In the spectra of the main shock and 
subregion (Fig. [5^ and b) the PAH emission bands at 

6.2 /im, 7.7 /jm and 8.6 fiui are far weaker compared to 
the 11.3 /Ltm bands, which when stron g are predominantl y 
produced by neutral PAH molecules ([Draine &: Lill200lh . 
In the shock subregion we find an upper limit flux for 
the 6.2 fim PAH of 0.47 x 10~ 17 Wm~ 2 and fluxes of 
0.87 x 10~ 17 Wm~ 2 and 2.49 x 10~ 17 Wm~ 2 for the 
7.7 and 11.3 //m bands respectively. This corresponds to 
a 7.7/zm/11.3/im PAH ratio of 0.35, very low compared 
to the median valu e found for the SINGS sample of 3.6 
([Smith et al.ll2007l ). 

The suppression of the 7.7/xm/11.3/im PAH ratio has 
been observed in AGN environ ments (e.g. lo w-luminosity 
AGN in the SINGS sample of lSmith I et all [2001 and is 
favored to be the result of selective destruction of PAH 
molecules small enough to emit at 7.7/xm. PAH process- 
ing in the shock due to larger molecules being less frag- 
ile than smaller ones is discussed further in Guillard et 
al. (2010). An alternative explanation is that the PAH 
molecules are chiefly large and neutral in the shock, pro- 
ducing enhanced 11.3/im emission in comparison to the 
smaller PAH molecules. A detailed comparison of the 
dust and PAH emission properties can be found in Guil- 
lard et al. (2010). 

Dust emission in the mid- and far-infrared can be 
used to infer the amount of star formation taking place 
(|Kennicutti 1 1 9981) . We measure a 24/j.m flux in this re- 
gion of 0.408 mjy, corresponding to a spectral luminosity 
[yh v ) of L24 = 1.38 x 10 Lq. This low luminosity is con- 
sistent with the weak MIR continuum (Fig. [SJd), arising 
from emission from VSGs heated by the UV radiation 
field. 

We can combine this with a measurement of the Ha 
emission in this region to obtain a star formation rate 
(SFR), given that they are complimentary (Ha trac- 
ing the young stellar population and 24/im as a mea- 
sure of dust-absorbed stellar light). We find Lhq = 

2.03 x 10 6 L Q , but caution that Ha emission in SQ is also 
the result of shock-excitation (|Xu et al.l l2003) and must 



be considered an upper limit for measuring star forma- 
ti on. When we combin e Lho with L24, using the relation 
of lCalzetti etHI ([2Q(51 . we find a SFR of <0.05 M Q yr" 1 . 

A further measure of star formation can be ob- 
tained from the P AH strength. Using the relation of 
| Houck et alJ (l2007f) . derived from the starburst sample of 
Bra ndl et al.l we can use the 7.7/im flux density 

as a measure of star formation. The shock subregion has 
a 7.7/xm flux density of 0.73 mjy which corresponds to a 
SFR of ~ 0.08 M Q yr" 1 , in good agreement with our pre- 
vious calculation, but also an upper limit as some PAH 
emission is contamination from known star formation re- 
gions in the group. A point of caution, however, is that 
this star formation indicator may be biased given that 
we detect suppressed 7.7/im emission compared to the 
11.3/im PAH. Comparing the Lpah(7 .7 /im) / L24 ratio to 
the galaxies in SINGS, we find t hat the value of ~ 0.17 is 
typical of star- forming galaxies (Ogle et al. 2009), which 
suggests that both measures of star formation are low, 
but self-consistent. 

The low upper limits for the SFR in the shock sug- 
gests that star formation is depressed in the shock apart 
from in the IGM starburst, SQ-A, which has a SFR of 
~ 1.2 5 M(n yr -1 at the velocity of the group (|Xu et al.l 
2003) . This would be consistent with a picture of molec- 
ular hydrogen being reheated by MHD shocks in the tur- 
bulent medium. 

Since the Jeans' mass increases with gas temperature 
and turbulence and shearing morions will prevent col- 
lapse, the cold molecular gas clouds may be too short- 
lived or undersized to f acilitate collapse and produce sig- 
nificant star formation (Guillard ct al. 2009]). The extent 
of a cold reservoir of molecular gas in the shock is a key 
consideration in this scenario, as discussed in Section lOl 

7. IMPLICATIONS FOR HIGH-Z GALAXY FORMATION 

The present observations suggest that molecular line 
cooling in dense clumps dissipates a significant fraction 
of the kinetic energy available in high-speed shocks. For 
reference, and to provide a glimpse of what this group 
might look like at high redshift when filling a single beam, 
we present in Figure [T3] the spectrum extracted for the 
entire group (i.e. including the three neighboring galax- 
ies and the shock). The total H2 line emission from the 
whole group exceeds 10 42 ergs _1 and is still the most 
dominant MIR coolant. Indeed the luminosity of the ro- 
tational H2 lines is sufficient that it could be detected 
at high redshift with future far-infrared or sub-mm in- 
strumentation like SPICA or SAFIR (see Appleton et al. 
2009). How likely is it, however, that high-speed shocks 
play a role in the assembly of galaxies? 

There is growing evi dence that galaxies at high 
redshift are turb ulent (|Forster Schreiber et al.l 120061 : 
IGenzel et al1[2008h and increasingly clumpy (e.g. Con- 
selice et al. 2005, Elmegreen & Elmegreen 2005). Indeed 
Bournaud & Elmegreen (2009) discuss the importance of 
the growth instabilities in massive gas clumps in form- 
ing disks at z > 1, and favor at least a large fraction 
of the clump systems being formed in smooth flows, per- 
haps similar to those discussed by Dekel et al. (2009). To 
what extent the build-up of these disks is truly "smooth" 
is not yet clear since the medium is likely to be a mul- 
tiphase one. In a more standard picture, gas flowing 
into the more massive dark halos will experience strong 
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Fig. 13. — Low resolution spectrum of the entire SQ region mapped by the IRS. This extraction is a rectangular box corre- 
sponding to an area of 94 kpc x 56 kpc centered on 22 h 35 m 59.97 s , +33°58'12.7", and includes the main shock, H 2 bridge, NGC 
7319, SQ-A and NGC 7318a and NGC 7318b. 



shocks, most likely in an inhomogenious medium (e. g. 
Greif et al. 2008) - thus it begs the question of how 
important H2 cooling may be in these different cases. 

Models of the collapse of the first structures predict 
that strong metal lines soon dominate the cooling over 
molecular hydrogen when the first stars pollute the en- 
vironment, ft is therefore usually assumed that, except 
at very early stages, molecular hydrogen is a minority 
coolant in gas that forms the first major structures (e.g. 
Bromm et al. 2001, Santoro and Shull 2006). However, 
our observations show that, under the right conditions, 
even in high metallicity environments, molecular hydro- 
gen can be extremely powerful - in this case dominating 
by a factor of ten over the usually powerful [Si n]34.82/im 
line in strong shocks. If there are situations at high red- 
shift where strong shocks propagate into a clumpy, multi- 
phase medium, then our observations imply that molecu- 
lar hydrogen cooling cannot be assumed to be negligible. 
On the other hand, this will not be a trivial problem. Our 
best model of Stephan's Quintet (|Guillard et alj l2009f) 
involves the formation of H2 in a complex multiphase 
turbulent medium in which shocks destroy dust in some 
places, but allow survival in others - thus encouraging 
H2 formation. 

In the early universe, this enhanced cooling, which has 
so far been neglected, will depend on the distribution 
and nature of the first dust grains, in concert with the 
formation, temperature and abundance of gas, and the 
feedback effects from the first stars and AGN. 



CONCLUSION 



In this paper we have presented the results of the mid- 
infrared spectral mapping of the Stephan's Quintet sys- 
tem using the Spitzer Space Telescope. We highlight here 
our five main conclusions: 

• The powerful H2 emission detected by 
lAppleton etaU (2006) surprisingly represents 
only a small fraction of the group-wide warm H2 
(with a lower limit luminosity of 1.1 x 10 42 ergs -1 
spread over ~ 975 kpc 2 ) that dominates the 
mid-infrared emission of the system. There is 
evidence for another shock-excited feature, the 
so-called H2 bridge between the main shock and 
NGC 7319, which is likely a remnant of past 
tidal interactions within the group. The spatial 
variation in the distribution of the H2 0-0 line 
ratios implies differences in temperature and 
excitation in the shocked system - this will be 
explored fully in Paper II. 



ratio in the main 



The global L(H2)/L x(0 .ooi-iokcv) 
shock is >= 3, and ^2.5 in the new "bridge" fea- 
ture. The results confirm that MIR H2 lines are 
a stronger coolant than X-ray emission over the 
shock structures, indicating a new cooling pathway 
seen on a large scale in SQ. This modifies the tra- 
ditional view that X-rays dominate cooling at all 
times in the later stages of compact group evolu- 
tion. Since H2 forms on the surfaces of dust grains, 
we expect dust emission associated with these re- 
gions, but a low intensity radiation field produces 
only weak emission at 24/im. 
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• Following earlier interpretations of nebular line ra- 
tios in the optical, we interpret infrared ionic lines 
within the framework of fast (V s > 100 kms" 1 ) 
ionizing shocks. Comparison between the [Neil], 
[Ne in] , [Si n] and [Fe n] line intensities implies that 
both silicon and iron are depleted onto dust. This 
result implies that dust is not destroyed in the 
shock. 

• Star formation in SQ is dominated by SQ-A and 
7318b-south, located at the extreme ends of the 
shock ridge seen at radio wavelengths, suggesting 
they are both shock triggered starbursts. However, 
regions dominated by warm H2 emission exhibit 
very low star formation rates, consistent with a 
turbulent model where H 2 is significantly reheated 
and cool clouds are too short-lived or undersized to 
collapse. 

• In SQ we observe the projected coexistence of 
[Sin] and H 2 being produced by ^200 km s -1 and 
^20kms _1 velocity shocks, respectively. Our ob- 
servational results are consistent with a model of 
a multiphase postsho ck medium produc ed by a 
galaxy- wide collision (jGuillard et al.ll2009| ). 
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APPENDIX A: THE SEYFERT II GALAXY, NGC 7319 

In this Section we discuss the results pertaining to NGC 7319, a Seyfert 2 galaxy (jDurretl fl994) lying to the East 
in the SQ group (see Fig. 1). The specific intensity contours for the H2 S(0) and H2 S(l) lines (Fig. 2a and b) show 
emission associated with the galaxy, seemingly connected to the rest of the group by the H2 "bridge" discussed in 
Section 13.11 Figure [3] shows that the n ucleus of NGC 7319 produces strong X-ray emission (jTrinchieri et all [2003) 
and is prominent at radio wavelengths (|Xu et al.ll2003f h However, we observe an offset between the peak of the H2 
emission near NGC 7319 and the Seyfert nucleus, which suggests that the "bridge" is a separate structure and not 
being excited by the AGN. 

The AGN in NGC 7319 does not have a well-collimat ed jet, but two extended lobes with compact hotspots, asy- 
metrically distributed along the minor axis of the galaxy fXanthopo ulos et al.ll2004 ). This structure runs NE/SW and 
its orientation compared to the H2 filament (which runs EW) is not consistent with causing the excitation of the H2 
bridge. This is also evident from the relatively weak power in the AGN, as inferred by the emission line diagnostics 
discussed below, and in the X-ray where it is only a factor of ^2 greater than the emission associated with the main 
shock and "bridge" . 

In Figure |4] we present the specific intensity contours of the [Feii]25.99/im (blended with [O iv]25.89/im), 
[S ni]33.48/im and [Si ii]34.82/im emission lines. Given the low spectral resolution of the SL and LL modules of 
Spitzer, we cannot distinguish between emission from [Fe 11] and [Oiv]. The AGN in NGC 7319 is li kely to produce 
both , with [Fe 11] emission likely originating in X-ray dissociation regions (XDRs) surrounding the AGN ( Mal onev et al.l 
1996). Prominent emission from [S iii]33.48/im and [Siii]34.82/mi are due to the high excitation conditions associated 
with the AGN. Dense PDRs and X-ray dominated regions, powered by AGN, show strong [Sin] emission at 34.82/mi, 
while [S ill] 33.48 [im emission acts as a tracer of Hn regions. 
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Fig. 14. — Rectangular extraction box for the Seyfert 2 galaxy NGC 7319 overlaid on the H2 0-0 S(l) line emission (a) and the 
corresponding low-resolution spectrum (b). The extraction box is centered on 22 h 36 m 03.62 s , +33°58'35.6". 

We now focus on the emission line properties of NGC 7319. Figure H4l shows the spectrum extracted from the galaxy 
and the measured properties of the line ratios are listed in Table [2] and [3l For the first time in the SQ group, we 
observe a spectrum which is no longer dominated by H2 emission, but instead the brightest lines arc the high-excitation 
[0iv]25.91/xm and [Fe n]25.98^m blended lines, as well as the [Si ii]34.81^m line. [Nev] is prominent at both 14.32^m 
and 24.30/im - a line typically seen in AGN. Also, unlike the majority of the extended shocked regions, there is a 
rising thermal continuum more typical of a starforming galaxy than a classical Seyfert galaxy, although as various 
studies have shown (Buchanan et al. 2006 ; Deo et al. 2007), Seyfert II galaxies exhibit a variety of MIR spectral 
characteristics at long wavelengths. NGC 7319's rising continuum is similar to that seen in the Seyfert II galaxy 
NGC 3079 (Deo et al. 2007), and likely represents a dominant starburst component in the Far-IR. Spitzer imaging 
reveals not only a bright nucleus, but also extended emission regions in the galaxy. 

Using the [Siii]34.81/rni/[Siii]33.48/xrn ratio as a probe of excitation sources, we find a ratio of ~1.85 which is low 
compared to the average value found for AGN galaxies (~ 2.9) in the Spitzer Infrared Nearby Galaxy Sample (SINGS), 
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but high compared to star-forming regions (~ 1.2) in the same sample (|Dale et al.l [20061). This suggests a relatively 
weak AGN. The [Nem] 15.56/im/[Neii] 12.81^m ratio is a measure of radiation field strength and the value of 0.97 
indicates a typical radiation field strength compared to other AGN (the sample of Weedman et al. 2005 shows a 
range of ~0.17 to 1 .9). The average electron dens ity, estimated from the [S in]18.71/jm/[S ml3 3.48^m ratio of ~ 0.56, 
is 100 - 200 cm -3 (|Martm-Hernandez et al.ll2002D , in the low-density limit for this diagnostic ([Smith et al.ll2009fl . 

The properties of the H 2 emission of NGC 7319 are limited to the two long wavelength lines observed by IRS-LL (the 
galaxy lies outside the region mapped by SL). As a result, the excitation diagram contains only two points allowing 
only an approximate idea of the H2 mass, since without the SL wavelength coverage to provide information on possible 
warmer components, we are likely to overestimate the temperature of the H2 by fitting a straight line to the points. 
Any warmer component would contribute to the flux of the 0-0 S(l) line thus leading to a reduction of the temperature 
(and increase in H2 mass) for any cooler component. However, to provide a guide, we estimate the temperature of H 2 
<=171±8 K (assuming the gas is in thermal equilibrium and thus an ortho-para ratio of 2.65) and a total H2 mass 
of 3.0±7x 10 8 Mq. We consider this a lower limit to the total warm mass if (as is likely) more than one component is 
present. Gao & Xu (2000) estimated the cold molecular hydrogen mass of NGC 7319 based on 12 CO (1-0) observations 
as > 3.6 x 10 9 Mp), a facto r of ro ughly ten greater than the warm H2 mass. Such a ratio is not atypical of large spiral 
galaxies iRigopoulou et al.l (|2002D . 

APPENDIX B: HIGH RESOLUTION SPECTRUM OF SHOCK CENTRE REANALYSED 

The high resolution spectrum at the center of the shock in SQ, obtained by lAppleton et al.l (|2006f ). has been 
reanalysed using the latest calibrations available for the IRS instrument (SSC pipeline version S17) and is included as 
Figure [T"5l This has shown that the H2 gas lies at a velocity of 6360(±100) km s _1 , between the velocity of the group 
(6600 kms -1 ) and the velocity of the intruder (5700 kms -1 ). This is consistent with a model of gas being accelerated 
by the shock, as well as the turbulence demonstrated by the broad linewidth of the H2 in the shock (860 kms -1 ). 
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Fig. 15.— Spitzer Short-High (SH) Spectrum, centred on 22 h 35 m 59.57 s , +33° 58' 1.8", reanalysed using the latest point source 
and extended source calibrations. 

The new SH spectrum shows that with the improved calibrations, the 11.3/im PAH feature is detected, although 
faint. This is consistent with what is found in the larger extractions that show enhanced 11.3/mi PAH emission 
compared to ionised PAH features emitting at 6.2/im and 7.7/xm. 
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APPENDIX C: X-RAY FLUXES IN EXTRACTED REGIONS 

In this section we present the X-ray fluxes in the extraction regions shown in Figure [7] The reanalysis of archival data 
is necessary to obtain accurate fluxes and luminosities for the X-ray emission over the vario us apertures matched to our 
spectral extractions. We use the XMM-Newton EPIC-pn data (see Trinchicri ct al. 2005, for observational details) to 
obtain the most sensitive measurements. A calibrated event file was generated and filtered using standard quality flags, 
and subsequently cleaned of background flares. A 0.3-2 keV image was then extracted and corrected for instrumental 
response. All point-like sources were masked out to 25" in the analysis of diffuse emission, and the local background 
level was estimated within a 5 x 9 arcmin 2 rectangular region away from the group core. For each region in Fig. 7, 
resulting background-subtracted photon count rates were converted to 0.3-2 keV and "bolometric" (0.001-10 keV) 
X-ray fluxes assuming an absorbed thermal plasma model of metallicity 0.4 Solar, with temperature as estimated from 
the map of O'Sullivan et al. (2009), and an absorbing Galactic Hi column density of Ah = 6.2 x 10 20 cm~ 2 . For 
NGC7319, which harbors an AGN, an absorbed power-law spectrum of photon index T = 1.7 was assumed instead. 
Results are listed in Table Q] 



TABLE 1 

X-ray Count Rates and Fluxes for Extracted Regions 



Region 


Count Rate (0.3-2 keV) 
(photons s - 1) 


Flux (0.3-2 keV) 
erg s~ 1 cm -2 


Bolometric Flux (0.001-10 keV) 
ergs -1 cm~ 2 


Adopted Spectral Model 


Main Shock 


0.102 


1.8 Xl0~ 13 


2.8 Xl0~ 13 


T = 0.7 keV 


Shock Sub-region 


0.017 


3.3 xl0~ 14 


5.1 Xl0~ 14 


T = 0.8 keV 


Bridge 


0.017 


3.0 xl0~ 14 


4.5 xl0~ 14 


T = 0.6 keV 


NGC 7319 


0.051 


1.1 xl0~ 13 


3.7 xl0~ 13 


r = 1.7 


SQ-A 


0.013 


2.4 xl0~ 14 


3.7 xl0~ 14 


T = 0.6 keV 



Note. — Statistical uncertainties on all flux measurements are < 10%. 



We note that for the main shock (with luminosity 1.9 xl0 41 ergs _1 in the 0-3-2 keV band) we are within a factor 
of 2 of the 0.5- 2 keV l uminosity ~3.1 x!0 41 er gs~ 1 obtained for a similar, but larger extraction of the shock by 
iTrinchieri et al.l (|2005| ). lO'Sullivan et al.l (|2Q09f > obtain a 0.5-2 keV surface brightness of 0.07 L© pc~ 2 in the main 
shock compared to our value of 0.1 L Q pc~ 2 for the 0.3-2 keV surface brightness. 
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TABLE 2 

Observed H2 Line Fluxes in units of 10 -17 W m -2 



Target Region 


Aperture 


H 2 0-0 S(0) 


H 2 0-0 S(l) 


H 2 0-0 S(2) 


H 2 0-0 S(3) 


H 2 0-0 S(4) 


H 2 0-0 S(5) 




(arcscc 2 ) 


A28.22/jm 


A17.03/jm 


A12.28jum 


A9.66^m 


A8.03/im 


A6.91/^m 


Main Shock 


2307 


3.09±0.19 


23.05±0.26 


9.10±0.38 


22.76±0.84 


2.5±1.0 


14.1±0.7 


Shock Sub-region 


242 


0.36±0.03 


3.13±0.02 


1.47±0.08 


3.60±0.08 


1.04±0.15 


2.37±0.30 


Bridge 


413 


0.72±0.04 


4.94±0.09 










NGC 7319 


1302 


2.06±0.26 


11.00±0.48 










SQ-A 


671 


1.83±0.05 


8.8±0.13 








3.7±0.93 



Note. — Upper limits arc the 3er values calculated from the RMS (root mean square) and the expected line profile width. 



TABLE 3 

Observed Fine Structure Line Fluxes in units of 10~ 17 W m~ 2 



Target Region 


[No 11] 
A12.81/im 


[Nc v] 
A14.32/jm 


[No in] 
A15.56/^m 


[Sin] 
A18.71^im 


[Nev] 
A24.32^m 


[Fcn] + [0 iv] 
A25. 99+25. 89pm 


[Sni] 
A33.48/im 


[Sin] 
A34.82/xm 


Main Shock 


5.31±0.30 


<0.9 


1.94±0.39 


0.76±0.13 


<0.14 


1.07±0.10 


1.87±0.15 


8.59±0.27 


Shock Sub-region 


1.04±0.09 


<0.18 


0.15±0.06 


<0.12 


<0.14 


0.15±0.03 


0.13±0.02 


1.23±0.04 


Bridge 






0.39±0.04 


<0.1 


<0.15 


0.29±0.04 


<0.22 


1.1±0.1 


NGC 7319 




2.89±0.64 


9.47±0.29 


2.81±0.40 


1.56±0.39 


10.9±0.27 


5.0±0.42 


9.24 ±0.49 


SQ-A 




<0.8 


1.29±0.16 


0.87±0.06 


<0.18 


<0.17 


1.54±0.10 


2.3±0.16 



Note. — Upper limits are the 3tr values calculated from the RMS and the expected line profile width. 



